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Preface 


This  report  presents  water  quality  monitoring  and  modeling  recommenda¬ 
tions  of  the  Lower  St.  Johns  River.  This  report  was  prepared  in  the  Environ¬ 
mental  Laboratory  (EL),  U.S.  Army  Engineer  Waterways  Experiment  Station 
(WES),  Vicksburg,  MS.  The  study  was  sponsored  by  the  U.S.  Army  Engi¬ 
neer  District,  Jacksonville,  and  was  funded  under  the  Military  Interdepartmen¬ 
tal  Purchase  Request  No.  RM-CW-92-0161,  dated  04  September  1992. 

The  Principal  Investigator  of  this  study  was  Mr.  Thomas  M.  Cole  of  the 
Water  Quality  and  Contaminant  Modeling  Branch  (WQCMB),  Environmental 
Processes  and  Effects  Division  (EPED),  EL.  This  report  was  prepared  by 
Mr.  Cole  under  the  direct  supervision  of  Dr.  Mark  Dortch,  Chief,  WQCMB, 
and  under  the  general  supervision  of  Mr.  Donald  L.  Robey,  Chief,  EPED, 
and  Dr.  John  W.  Keeley,  Director,  EL.  Technical  reviews  by  Dr.  Dortch  and 
Ms.  Dorothy  H.  Tillman  are  gratefully  acknowledged. 

At  the  time  of  publication  of  this  report.  Director  of  WES  was 
Dr.  Robert  W.  Whalin.  Commander  of  WES  was  COL  Bruce  K.  Howard, 
EN. 

This  report  should  be  cited  as  follows; 

Cole,  T.  M.  (1995).  "Review  of  water  quality  monitoring 
and  recommendations  for  water  quality  modeling  of  the  Lower 
St.  Johns  River,"  Miscellaneous  Paper  EL-95-3,  U.S.  Army 
Engineer  Waterways  Experiment  Station,  Vicksburg,  MS. 


The  contents  of  this  report  are  not  to  be  used  for  advertising,  publication, 
or  promotional  purposes.  Citation  of  trade  names  does  not  constitute  an 
official  endorsement  or  approval  of  the  use  of  such  commercial  products. 


Introduction 


The  St.  Johns  River  Water  Management  District  (SJRWMD)  has  plans  to  devel¬ 
op  a  comprehensive  modeling  system  for  the  Lower  St.  Johns  River  Estuary 
(LSJRE).  The  modeling  system  will  be  used  to  provide  information  for  issues 
concerning  watershed  management,  water  supply,  pollutant  loadings,  and  eutro¬ 
phication.  Thus,  the  modeling  system  will  require  models  of  the  watershed, 
hydrodynamics,  sediment  transport  processes,  and  water  quality  including 
algal/nutrient/dissolved  oxygen  interactions  and  contaminants.  The  water  quali¬ 
ty  model  is  an  important  component  of  the  system.  In  light  of  this,  the 
SJWRMD  has  requested  a  review  of  the  present  water  quality  monitoring  pro¬ 
gram  to  insure  that  proper  and  sufficient  data  are  being  collected  to  support 
future  water  quality  modeling  efforts. 

This  report  reviews  present  water  quality  data  gathering  efforts,  recommends 
additional  data  efforts,  and  recommends  and  discusses  a  water  quality  model 
suitable  to  meet  the  objectives  of  the  SJRWMD. 
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Water  Quality  Monitoring 


Data  provided  by  SJRWMD  for  review  were  limited  in  scope.  These  data 
included  a  list  of  parameters  presently  sampled  by  the  SJRWMD,  sampling 
station  locations  used  by  SJRWMD  and  other  agencies,  values  of  typical  water 
quality  parameters  plotted  spatially  from  just  downstream  of  the  Ocklawaha 
River  to  the  estuary  mouth  for  one  sampling  event,  and  average  values  from 
1983  to  1992  for  three  sections  of  the  St.  Johns  River.  Therefore,  a  compre¬ 
hensive  review  of  existing  data  is  not  possible  in  this  report.  The  following 
provides  a  limited  discussion  of  the  provided  water  quality  data  and  discusses  in 
general  terms  water  quality  data  requirements  necessary  to  support  a  water 
quality  model  of  the  LSJR. 


Existing  Monitoring 

Sampled  Parameters 

Table  1  lists  the  parameters  routinely  sampled  for  by  SJRWMD.  In  addition, 
metals  including  calcium,  magnesium,  chromium,  zinc,  lead,  cadmium,  and 
copper  are  also  monitored. 

Spatial  Monitoring 

The  SJRWMD  has  an  extensive  water  quality  sampling  program  in  place. 
Information  provided  by  SJRWMD  shows  that  approximately  40  stations  are 
routinely  sampled  throughout  the  Lower  St.  Johns  River  by  the  SJRWMD. 
Approximately  50  additional  stations  are  sampled  by  other  agencies.  No  infor¬ 
mation  was  supplied  as  to  what  water  quality  parameters  were  sampled  at  the 
additional  stations.  Major  tributaries  appear  to  be  adequately  represented  allow¬ 
ing  for  development  of  river  loading  estimates.  No  information  was  provided 
on  the  vertical  extent  of  sampling. 

Temporal  Monitoring 

Water  quality  parameters  are  typically  sampled  once  a  month.  In  addition, 
synoptic  sampling  events  the  length  of  the  lower  river  are  conducted  approxi¬ 
mately  twice  a  year.  No  information  was  provided  as  to  what  parameters  are 
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sampled  for  during  the  synoptic  events.  It  is  assumed  that  the  parameters 
include  those  listed  in  Table  1. 


Table  1 

Sampled  Water  Quality  Parameters 

Parameter 

Frequency 

Temperature 

monthly 

pH 

monthly 

Dissolved  oxygen 

monthly 

Ammonium-N 

monthly 

Nitrate +  Nitrite-N 

monthly 

Total  Kjeldahl  N 

monthly 

Total  phosphorus 

monthly 

Orthophosphate 

monthly 

Dissolved  silica 

monthly 

Total  suspended  solids 

monthly 

Total  dissolved  solids 

monthly 

Sulfate 

monthly 

Chloride 

monthly 

Chlorophyll  A,BrC,  and  phaephyton 

monthly 

Total  organic  carbon 

monthly 

Alkalinity 

monthly 

Turbidity 

monthly 

Color 

monthly 

Determination  of  Loadings 

A  key  to  any  successful!  water  quality  modeling  effort  is  an  accurate  determina¬ 
tion  of  the  loadings  to  the  system.  Without  the  correct  loadings,  calibration  of 
the  model  is  difficult  and  results  of  management  scenarios  will  be  suspect  at 
best.  It  is  therefore  crucial  to  obtain  accurate  loading  estimates.  Stations  at 
tributary  mouths  are  presently  sampled  monthly.  This  is  not  sufficient  to  deter¬ 
mine  accurate  estimates  of  tributary  loadings.  A  sampling  procedure  for  esti¬ 
mating  loads  from  tributary  loads  will  be  discussed  in  more  detail  in  the  next 
section. 
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Recommended  Water  Quality  Monitoring 


Table  2  lists  the  recommended  parameters  necessary  for  supporting  a  state-of- 
the-art  water  quality  model  of  the  Lower  St.  Johns  River.  Of  primary  impor¬ 
tance  are  particulate  (POC)  and  dissolved  organic  carbon  (DOC),  primary 
production,  and  sediment  oxygen  demand  and  nutrient  fluxes.  Dissolved  organ¬ 
ic  carbon  plays  a  key  role  in  eutrophication  related  issues  in  the  LSJR.  Based 
on  color  measurements  and  visual  inspection  of  the  LSJR,  it  appears  that  pro¬ 
ductivity  in  the  fresh  water  portion  of  the  river  is  limited  by  light  penetration 
due  primarily  to  DOC.  Discussions  with  personnel  from  SJRWMD  suggest  that 
loadings  of  DOC  are  primarily  from  major  tributaries  and  wetland  areas  in  the 
upper  portion  of  the  river.  Therefore,  it  is  imperative  that  accurate  regressions 
of  DOC  with  flow  be  developed  for  estimating  DOC  loadings  to  the  system. 
Total  inorganic  carbon,  along  with  pH  and  alkalinity,  will  allow  modeling  of 
the  pH-carbonate  system. 

Also  of  importance  are  determinations  of  phytoplankton  productivity  and  respi¬ 
ration  rates  and  bioassays  for  nutrient  limitations.  Water  quality  models  are 
capable  of  computing  nutrient  and  light  limitations  and  bioassays  would  be 
helpful  in  comparing  model  predictions.  Productivity  and  respiration  rates  are 
necessary  for  accurate  characterization  of  carbon  assimilation  by  phytoplankton. 

It  is  also  critical  to  obtain  accurate  estimates  of  sediment  oxygen  demand  and 
nutrient  fluxes  for  proper  calibration  of  a  sediment  model  or  for  proper  specifi¬ 
cation  of  boundary  values  if  a  sediment  model  is  not  used.  In  addition,  proper 
characterization  of  the  sediments  including  particulate  organic  carbon,  phospho¬ 
rus,  nitrogen,  and  inorganic  phosphorus  is  necessary  for  sediment  model  cali¬ 
bration. 

Total  carbon,  nitrogen,  phosphorus,  and  silica  are  necessary  to  insure  that 
loadings  during  model  calibration  are  correctly  accounted  for.  Care  should  be 
taken  during  analysis  so  that  no  single  component  is  greater  than  the  total 
nutrient.  If  this  occurs,  samples  should  be  immediately  analyzed  again. 

Monthly  sampling  frequency  for  all  in-pool  parameters  is  generally  sufficient  in 
the  early  stages  of  a  sampling  program.  However,  as  knowledge  of  the  system 
behavior  increases,  sampling  frequency  during  the  most  dynamic  periods  may 
have  to  be  increased  to  better  capture  the  dynamics  of  the  system.  For  exam¬ 
ple,  it  may  be  necessary  to  sample  more  frequently  during  shifts  in  phytoplank¬ 
ton  assemblages  in  order  to  gain  a  better  understanding  of  phytoplankton  suc¬ 
cession.  It  is  critical  to  have  an  understanding  of  the  system  in  order  to  accu¬ 
rately  represent  the  processes  during  a  modeling  effort.  In  addition,  most 
modeling  studies  give  direction  to  future  sampling  efforts  by  pinpointing  def- 
iciences  in  the  existing  data. 

Concerning  contaminant  modeling,  it  is  hard  to  make  specific  recommendations 
without  knowing  the  contaminant  of  concern.  However,  the  following  parame- 
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Table  2 

Recommended  Monitoring  Parameters 

Water  Column 

Temperature 

Dissolved  phosphate 

Salinity 

* 

Particulate  inorganic  phosphorus 

Light  extinction 

♦ 

Dissolved  organic  phosphorus 

Dissolved  oxygen 

* 

Particulate  organic  phosphorus 

Ammonium 

* 

Total  phosphorus 

Nitratenitrlte 

Dissolved  silica 

* 

Dissolved  organic  nitrogen 

* 

Particulate  biogenic  silica 

♦ 

Particulate  organic  nitrogen 

* 

Total  silica 

* 

Total  Nitrogen 

Chlorophyll 

pH 

* 

Total  inorganic  carbon 

Alkalinity 

Sediments 

* 

Particulate  organic  carbon 

* 

Particulate  inorganic  phosphorus 

* 

Particulate  organic  nitrogen 

* 

Particulate  organic  phosphorus 

Processes/Fluxes 

* 

Primary  production 

* 

Sediment  nitrate  flux 

Respiration 

* 

Sediment  phosphorus  flux 

* 

Sediment  oxygen  demand 

«• 

Sediment  silica  flux 

* 

Sediment  ammonium  flux 

Sediment  COD  flux* 

♦ 

Nutrient  limitation  bioassays 

■N- 

recommended  additional  sampling  and  analyses 

ters  should  be  measured  for  any  contaminant  modeling  study.  In  the  sediments, 
bulk  and  dry  density,  grain  size,  TOC  and  DOC,  longterm  burial  velocities, 
and  contaminant  concentrations  should  be  measured.  Longterm  burial  rates  are 
typically  determined  from  strontium  core  datings.  In  the  water  column,  total 
and  dissolved  contaminant  concentrations,  TOC  and  DOC,  and  settling  velocity 
should  be  measured.  In  addition,  sorption  and  degradation  values  should  be 
determined  from  either  the  literature  or  from  laboratory  experiments. 
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Recommended  Water  Quality  Model 

The  recognized  state-of-the-art  water  quality  model  with  respect  to  eutrophica¬ 
tion  related  issues  is  CE-QUAL-ICM  (Cerco  and  Cole,  1993).  The  model  is 
similar  in  concept  to  WASP  (Ambrose,  et.al.,  1986)  and  was  developed  over 
the  last  six  years  by  the  Water  Quality  and  Contaminant  Modeling  Branch  at 
Waterways  Experiment  Station  for  application  to  Chesapeake  Bay.  The  model 
is  presently  being  used  by  the  U.S.  Army  Engineer  Waterways  Experiment 
Station  (WES)  to  evaluate  the  effectiveness  of  nutrient  reduction  strategies  on 
eutrophication  in  the  Bay,  and  has  also  been  used  to  address  management  issues 
in  the  New  York  Bight,  Indian  River-Reheboth  Bay,  Delaware,  LA  Long 
Beach,  and  Lower  Green  Bay.  CE-QUAL-ICM  is  the  first  3-D  water  quality 
model  shown  to  be  capable  of  addressing  water  quality  issues  spanning  multiple 
decades. 

Capabilities 

Mass  Balance.  Water  quality  modeling  can  be  viewed  simply  as  mass-balance 
modeling  wherein  the  model  keeps  track  of  constituent  mass  over  time  as  it  is 
affected  by  external  boundary  conditions  and  internal  kinetic  transformations. 
CE-QUAL-ICM  has  been  shown  to  account  for  constituent  mass  within  the 
accuracy  of  a  given  computer.  This  is  the  most  fundamental  requirement  of 
any  water  quality  model  -  the  water  quality  transformations  and  numerical 
solution  scheme  must  maintain  mass  balances. 

Hydrodynamic  Linkage.  CE-QUAL-ICM  is  a  finite-volume  model.  Theoreti¬ 
cally,  the  model  can  link  to  hydrodynamic  models  of  any  dimension  or  combi¬ 
nation  of  dimensions  using  either  structured  or  unstructured  grids.  However, 
certain  numerical  solution  techniques  in  the  model  complicate  the  linkage  to 
unstructured  grid,  hydrodynamic  models.  This  will  be  discussed  in  more  detail 
in  the  Limitations  section. 

Transport.  CE-QUAL-ICM  is  a  stand  alone  water  quality  code  decoupled 
from  the  hydrodynamics.  Output,  specifically  flows  and  diffusivities  across  cell 
faces,  from  a  chosen  hydrodynamic  model  is  used  to  drive  the  transport  pro¬ 
cesses.  Diffusivities  can  be  optionally  specified  in  the  water  quality  model. 
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The  model  is  capable  of  addressing  long-term  water  quality  issues  in  one,  two, 
or  three  dimensions,  or  in  any  combination.  Thus,  a  variety  of  hydrodynamic 
models  can  be  used  to  provide  transport  information  required  by  CE-QUAL- 
ICM. 

A  major  advancement  incorporated  in  the  model  is  a  subcomponent  which 
calculates  residual  currents  from  the  hydrodynamic  model  averaged  over  a 
given  time  period  (typically  a  tidal  cycle).  This  allows  for  longer  timesteps  in 
the  water  quality  model  while  still  maintaining  the  transport  properties  of  the 
hydrodynamic  model  and  contributes  to  the  model’s  ability  to  perform  longterm 
simulations.  The  calculation  of  residual  currents  is  presently  incorporated  into 
the  hydrodynamic  model  CH3D  (Dortch,  1992).  Additional  effort  would  be 
required  to  provide  residual  flows  from  a  different  hydrodynamic  model. 

Another  important  advancement  incorporated  into  CE-QUAL-ICM  is  a  higher- 
order  numerical  transport  scheme,  QUICKEST  (Leonard,  1979),  which  greatly 
reduces  numerical  diffusion  during  constituent  transport.  It  is  well  known  that 
upwind  differencing  produces  numerical  dispersion  which  can  often  be  much 
greater  than  the  physical  dispersion.  Because  of  numerical  dispersion,  upwind 
differencing  cannot  resolve  sharp  gradients  which  are  typically  encountered 
when  modeling  water  quality.  The  transport  solution  scheme  in  CE-QUAL- 
ICM  can  use  either  upwind  differencing  or  QUICKEST.  The  preferred  method 
is  QUICKEST.  The  Chesapeake  Bay  Study  showed  that  observed  vertical 
gradients  in  salinity  could  not  be  reproduced  in  the  hydrodynamic  model  with¬ 
out  the  inclusion  of  a  higher-order  transport  scheme  (Johnson,  1991). 

Vertical  transport  is  differentiated  from  horizontal  transport  in  the  model. 
Vertical  diffusion  is  always  handled  implicitly  in  time  which  removes  vertical 
diffusion  from  numerical  stability  criteria.  Vertical  advection  employs  a  time- 
weighted  scheme  which  can  vary  between  fully  explicit,  fiilly  implicit,  or  a 
combination  of  the  two.  The  fully  implicit  option  also  eliminates  vertical  ad¬ 
vection  from  the  transport  stability  criteria.  The  combination  of  fully  implicit 
vertical  advection  and  diffusion  contributes  to  the  model’s  ability  to  conduct 
longterm  simulations  (ie.,  decades)  required  when  evaluating  nutrient  reduction 
scenarios  or  contaminant  modeling. 

The  model  also  includes  an  algorithm  which  adjusts  the  timestep  based  on 
estimates  of  numerical  stability  criteria  required  by  the  numerical  solution  of  the 
advection-diffusion  equation.  Thus,  the  timestep  automatically  decreases  during 
periods  of  high  constituent  fluxes  and  increases  during  low  periods.  Previous¬ 
ly,  the  timestep  in  most  water  quality  models  was  set  to  the  lowest  value  en¬ 
countered  during  a  simulation  to  maintain  numerical  stability.  Automatic  model 
adjustment  of  the  timestep  is  referred  to  as  autostepping  and  has  a  dramatic 
effect  on  reducing  the  computational  requirements  of  the  model. 

Eutrophication  Interactions.  The  model  contains  a  full  suite  of  water  quality 
constituents  used  to  address  eutrophication  related  issues.  A  complete  descrip¬ 
tion  of  the  state  variables  and  their  interactions  is  given  in  Appendix  A.  This 
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appendix  is  provided  primarily  to  allow  the  reader  to  determine  if  the  important 
processes  and  proper  formulations  for  the  LSJR  are  included  in  the  model. 
This  appendix  was  taken  from  a  draft  report  of  the  Chesapeake  Bay  study 
(Cerco  and  Cole,  1993). 

Sediment  Submodel.  CE-QUAL-ICM  contains  a  predictive  sediment  submodel 
which  represents  the  accumulation  and  decay  or  diagenisis  of  organic  matter  in 
the  sediments  and  its  effect  on  water  column  oxygen  demand  and  nutrient  fluxes 
from  the  sediments  (DiToro  and  Fitzpatrick,  1993).  The  present  sediment 
formulation  is  applicable  mainly  to  saline  conditions  in  the  water  column.  A 
freshwater  formulation  which  includes  iron  and  manganese  will  be  incorporated 
into  the  model  in  the  future. 

Limitations.  It  was  previously  stated  that  CE-QUAL-ICM  is  the  state-of-the- 
art  in  water  quality  modeling.  However,  it  does  not  contain  all  of  the  state-of- 
the-art  capabilities  available  today  in  the  field  of  water  quality  modeling.  It 
should  be  noted  that  no  other  water  quality  model  contains  complete  state-of- 
the-art  capabilities  either.  The  model  was  initially  developed  to  address  specific 
management  issues  for  the  Chesapeake  Bay  Study,  but  the  coding  of  the  model 
is  such  that  additional  capabilities  can  be  incorporated  as  needed.  In  addition, 
the  algorithms  in  the  code  are  applicable  but  not  specific  to  Chesapeake  Bay. 

Experience  has  shown  that  nearly  every  new  application  of  any  water  quality 
model  requires  inclusion  of  additional  processes  or  reformulation  of  existing 
processes  to  adequately  reproduce  prototype  behavior.  CE-QUAL-ICM  will 
evolve  over  time  as  new  applications  require  additional  capabilities  or  reformu¬ 
lations.  The  most  important  limitations  of  the  model  are  discussed  below. 

Hydrodynamic  Linkage.  The  only  hydrodynamic  model  used  to  drive  the 
transport  processes  in  previous  CE-QUAL-ICM  applications  is  CH3D-WES 
(Johnson,  1991).  CH3D-WES  is  discussed  in  the  companion  report  and  will  be 
discussed  only  briefly  in  this  report.  The  discretization  grid  is  boundary  fitted, 
curvilinear,  structured,  and  space  staggered.  The  vertical  grid  can  be  either 
sigma-stretched  where  the  number  of  layers  in  the  grid  are  fixed  and  layer 
thicknesses  vary  temporally  throughout  the  water  column,  or  a  Z-grid  represen¬ 
tation  where  only  the  surface  layer  varies  in  thickness  over  time.  CE-QUAL- 
ICM  is  capable  of  linking  to  either  vertical  grid.  Most  importantly,  the  higher- 
order  transport  scheme  in  CE-QUAL-ICM  requires  a  structured  grid  representa¬ 
tion.  CH3D  meets  this  requirement. 

Theoretically,  finite-element  hydrodynamic  models  using  unstructured  grids  can 
be  used  to  drive  the  transport  processes  in  CE-QUAL-ICM.  However,  finite- 
element,  unstructured  grids  have  not  previously  been  used  with  CE-QUAL- 
ICM.  Plans  have  been  made  to  incorporate  the  ability  to  link  CE-QUAL-ICM 
to  unstructured  grids,  but  it  must  be  emphasized  that  research  would  be  involv¬ 
ed  in  linking  the  water  quality  model  to  a  finite-element,  hydrodynamic  model. 
The  following  considerations  are  involved  with  linking  CE-QUAL-ICM  to  a 
finite-element  hydrodynamic  model. 
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Finite-element  hydrodynamic  models  conserve  mass  globally,  but  do  not  con¬ 
serve  mass  locally.  As  stated  previously,  the  fundamental  principle  in  water 
quality  modeling  involves  keeping  track  of  constituent  mass  through  external 
boundary  interactions  and  internal  kinetic  transformations.  Presently,  it  is  not 
known  what  affect  the  inability  of  a  hydrodynamic  model  to  locally  conserve 
mass  will  have  on  the  ability  of  CE-QUAL-ICM  to  reproduce  behavior  in  the 
prototype.  It  is  believed  that  finite-element  hydrodynamic  models  can  be  cor¬ 
rected  to  conserve  mass  locally,  but,  as  stated  previously,  research  would  be 
required  to  establish  this. 

More  importantly,  the  present  higher-order  transport  scheme  (QUICKEST)  in 
CE-QUAL-ICM  requires  three  spatial  points  to  interpolate  constituent  fluxes 
across  cell  faces.  Linking  CE-QUAL-ICM  to  an  unstructured  hydrodynamic 
model  would  require  the  development  of  a  two-point,  higher-order  transport 
scheme  since  the  three  points  required  by  the  present  scheme  would  not  be 
readily  discernable.  Such  schemes  are  available,  but  again,  research  would  be 
involved  in  implementing  the  two-point,  higher-order  transport  scheme.  Addi¬ 
tionally,  the  new  transport  scheme  must  be  shown  to  reproduce  the  transportive 
properties  of  the  hydrodynamic  solution  scheme. 

Future  Capabilities 

Numerous  improvements  to  CE-QUAL-ICM  are  planned  in  the  future.  They 
are  discussed  in  the  following  sections. 

Sediment  Submodel  w/  Iron  and  Manganese  Kinetics.  The  current  sediment 
model  was  developed  to  compute  carbon  diagenesis  and  sediment-water  column 
fluxes  for  Chesapeake  Bay  which  is  predominately  saltwater.  Development  is 
underway  to  make  the  sediment  submodel  applicable  to  freshwater  sediments. 
The  main  differences  between  the  freshwater  and  saltwater  formulations  are  that 
sulfides  are  the  primary  source  of  oxygen  demand  in  saltwater  and  methane, 
iron,  and  manganese  along  with  sulfides  are.  the  primary  sources  of  sediment 
oxygen  demand  in  freshwater. 

Carbonate-pH  System.  A  carbonate-pH  system  with  non-conservative  alkalini¬ 
ty  is  being  incorporated  into  the  model.  A  fully  tested  system  is  expected  by 
1995. 

Zooplankton.  The  initial  water  quality  formulations  in  CE-QUAL-ICM  explic¬ 
itly  included  the  effects  of  zooplankton  by  their  inclusion  as  state  variables. 
During  model  calibration,  zooplankton  were  removed  from  the  formulations  and 
their  effects  were  included  in  the  algal  mortality  and  nutrient  recycling  terms. 
Plans  are  to  reintroduce  zooplankton  into  the  model  as  state  variables  for  Phase 
2  of  the  Chesapeake  Bay  Study. 
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Macrophytes.  Macrophyte/nutrient  interactions  are  not  presently  included  in 
CE-QUAL-ICM.  Phase  2  of  the  Chesapeake  Bay  Study  will  incorporate  macro¬ 
phytes  as  a  state  variable.  This  is  presently  an  active  area  of  research  in  the 
modeling  community  and  improvements  will  be  incorporated  in  the  future  as  a 
better  understanding  of  the  processes  are  gained. 

Contaminants.  The  contaminant  components  of  the  EPA’s  WASP  model 
(Ambrose,  et.  al.,  1986)  are  presently  being  incorporated  into  the  model.  No 
contaminant  studies  have  been  conducted  with  CE-QUAL-ICM;  however,  a 
water  quality  and  contaminant  modeling  study  of  San  Diego  Bay  is  scheduled  to 
start  in  1994. 

Sediment  Transport.  Sediment  transport  is  presently  not  in  the  model.  It 
should  be  relatively  easy  to  incorporate  sediment  transport  given  sediment 
transport  output  from  an  appropriate  hydrodynamic  model.  Again,  sediment 
transport  will  probably  be  incorporated  into  the  model  during  the  San  Diego 
Bay  study. 

Pre-  and  Postprocessors.  One  of  the  most  serious  shortcomings  of  CE-QUAL- 
ICM  is  a  lack  of  integrated  pre-  and  postprocessors.  Presently,  the  develop¬ 
ment  of  input  files  requires  a  hodgepodge  of  FORTRAN  programs  whose 
workings  are  largely  undocumented.  The  same  is  true  for  postprocessing  of 
output.  Graphical  display  of  model  output  consists  of  a  series  of  FORTRAN 
programs  which  use  high  level  graphics  subroutine  calls  from  a  commercial 
package  called  DISSPLA.  There  are  plans  to  develop  integrated  processors 
based  on  a  graphical  user  interface.  However,  given  the  current  level  of  fund¬ 
ing,  development  will  take  several  years  to  complete. 

Computer  Requirements 

Hardware.  CE-QUAL-ICM  has  previously  been  run  on  80386/80486  PC’s, 
RISC  workstations,  and  a  Cray  Y-MP  and  C-90.  All  of  the  code  developement 
for  Chesapeake  Bay  was  performed  on  a  25  Mhz  80386  PC  while  model  cali¬ 
bration  and  scenario  runs  were  performed  on  a  Cray  Y-MP.  Given  the  rapid 
advancements  in  the  state-of-the-art  in  hardware  for  computationally  intensive 
floating-point  code,  it  is  hard  to  estimate  future  hardware  requirements  for 
management  scenarios.  It  is  almost  a  certainty  that  the  model  will  run  in  a 
reasonable  time  on  workstations  costing  less  than  lOK  in  five  years. 

Software.  Present  software  requirements  are  a  FORTRAN  compiler  and  a 
commercial  software  package  called  DISSPLA  for  processing  of  model  output. 
For  PC’s,  the  compiler  must  be  capable  of  addressing  32-bits  in  DOS’s  protect¬ 
ed  mode.  The  trend  in  workstations  is  toward  64-bit  CPU’s  which  deliver  the 
same  precision  as  Cray  supercomputers.  Machines  with  64-bits  of  precision  are 
not  a  requirement  of  the  model,  but  they  increase  the  accuracy  of  the  computa¬ 
tions. 
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Appendix  A 
CE-QUAL-ICM  Kinetics 


The  central  issues  in  the  water  quality  model  are  primary  production  of  carbon 
by  algae  and  dissolved  oxygen  concentrations.  Primary  production  provides  the 
energy  required  by  the  ecosystem  to  function.  Excessive  primary  production  is 
detrimental  since  its  decomposition  in  the  water  and  sediments  consumes  oxy¬ 
gen.  Dissolved  oxygen  is  necessary  to  support  the  life  functions  of  higher 
organisms  and  is  considered  an  indicator  of  the  "health”  of  estuarine  systems. 
In  order  to  predict  primary  production  and  dissolved  oxygen,  a  large  suite  of 
model  state  variables  is  necessary  (Table  A-1), 


Table  A-1 

Water  Quality  Model  State  Variables 

temperature 

Salinity 

total  active  metal 

cyanobacteria 

diatoms 

green  algae 

dissolved  organic  carbon 

labile  particulate  organic  carbon 

refractory  particulate  organic  carbon 

ammonium 

nitrate 

dissolved  organic  nitrogen 

labile  particulate  organic  nitrogen 

refractory  particulate  organic  nitrogen 

total  phosphate 

dissolved  organic  phosphorus 

labile  particulate  organic  phosphorus 

refractory  particulate  organic  phos¬ 
phorus 

chemical  oxygen  demand 

dissolved  oxygen 

particulate  biogenic  silica 

available  silica 
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Conservation  Of  Mass  Equation 


The  foundation  of  CE-QUAL-ICM  is  the  solution  to  the  three-dimensional 
mass-conservation  equation  for  a  control  volume.  The  control-volume  approach 
was  selected  to  allow  maximum  flexibility  in  linking  CE-QUAL-ICM  to  alter¬ 
nate  hydrodynamic  models.  For  each  state  variable,  the  conservation  of  mass 
is: 


6V,Ci 

“5t~ 


H 


j=i 


6x 


(A-1) 


where 

Vi  =  volume  (m^) 

Ci  =  concentration  (g  m'^) 

Qj  =  volumetric  flow  across  flow  face  j  (m^  sec  ') 

C'j  =  concentration  in  flow  across  flow  face  j  (g  m'^) 

Aj  =  area  of  flow  face  j  (m^) 

Dj  =  diffusion  coefficient  at  flow  face  j  (m^  sec  ’) 
n  =  number  of  flow  faces  attached  to  control  volume 
Sj  =  control  volume  external  loads  and  kinetic  sources/sinks  (g  sec  ') 
t,  X  =  temporal  and  spatial  coordinates 


Solution  to  the  mass-conservation  equation  uses  the  QUICKEST  finite  differ¬ 
ence  algorithm  (Leonard,  1979)  for  horizontal  advection  and  diffusion  and  a 
Crank-Nicolson  scheme  for  vertical  advection  and  diffusion. 


State  Variables 

Algae 

Algae  are  grouped  into  three  model  classes:  cyanobacteria,  diatoms,  and  greens. 
The  grouping  is  based  upon  the  distinctive  characteristics  of  each  class  and 
upon  the  significant  role  they  play  in  the  ecosystem.  Cyanobacteria,  commonly 
called  blue-green  algae,  are  characterized  by  their  abundance  (as  picoplankton) 
in  saline  water  and  by  their  bloom-forming  characteristics  in  fresh  water. 
Cyanobacteria  are  unique  in  that  some  species  fix  atmospheric  nitrogen.  The 
cyanobacteria  in  the  model  are  bloom-forming  species  that  do  not  fix  nitrogen. 
They  are  characterized  as  having  negligible  settling  velocity  and  subject  to  low 
predation  pressure.  Diatoms  are  distinguished  by  their  requirement  of  silica  as 
a  nutrient  to  form  cell  walls  with  high  settling  velocities.  Settling  of  spring 
diatom  blooms  to  the  sediments  may  be  a  significant  source  of  carbon  for  sedi- 
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ment  oxygen  demand.  Algae  that  do  not  fall  into  the  preceding  two  groups  are 
lumped  into  the  heading  of  green  algae.  Green  algae  settle  at  a  rate  interme¬ 
diate  between  cyanobacteria  and  diatoms  and  are  subject  to  greater  grazing 
pressure  than  cyanobacteria. 

Organic  Carbon 

Three  organic  carbon  state  variables  are  included:  dissolved,  labile  particulate, 
and  refractory  particulate.  Labile  and  refractory  distinctions  are  based  upon  the 
time  scale  of  decomposition.  Labile  organic  carbon  decomposes  on  a  time  scale 
of  days  to  weeks  while  refractory  organic  carbon  requires  more  time.  Labile 
organic  carbon  decomposes  rapidly  in  the  water  column  or  the  sediments. 
Refractory  organic  carbon  decomposes  slowly,  primarily  in  the  sediments,  and 
contributes  to  sediment  oxygen  demand  years  after  deposition. 

Nitrogen 

Nitrogen  is  divided  into  organic  and  mineral  fractions.  Organic  nitrogen  state 
variables  are:  dissolved  organic  nitrogen,  labile  particulate  organic  nitrogen, 
and  refractory  particulate  organic  nitrogen.  Two  mineral  nitrogen  forms  are 
modeled  -  ammonium  and  nitrate/nitrite.  Both  forms  are  available  for  algal 
nutrient  requirements  although  ammonium  is  preferred  from  thermodynamic 
considerations.  The  primary  reason  for  distinguishing  the  two  is  that  ammoni¬ 
um  is  oxidized  by  nitrifying  bacteria  into  nitrate.  This  oxidation  can  be  a 
significant  sink  of  oxygen  in  the  water  column  and  sediments.  An  intermediate 
in  the  complete  oxidation  of  ammonium,  nitrite,  also  exists.  Nitrite  concentra¬ 
tions  are  usually  much  less  than  nitrate  and,  for  modeling  purposes,  are  com¬ 
bined  with  nitrate.  The  nitrate  state  variable  represents  the  sum  of  nitrate  plus 
nitrite. 

Phosphorus 

As  with  carbon  and  nitrogen,  organic  phosphorus  is  considered  in  three  states: 
dissolved,  labile  particulate,  and  refractory  particulate.  Only  a  single  mineral 
form,  total  phosphate,  is  considered.  Total  phosphate  exists  in  several  states 
within  the  model:  dissolved  phosphate,  phosphate  sorbed  to  inorganic  solids, 
and  phosphate  incorporated  in  algal  cells.  Equilibrium  partition  coefficients  are 
used  to  distribute  the  total  among  the  three  states. 
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Silica 


Silica  is  divided  into  two  state  variables:  available  silica  and  particulate  biogenic 
silica.  Available  silica  is  primarily  dissolved  and  can  be  used  by  diatoms. 
Particulate  biogenic  silica  cannot  be  used.  In  the  model,  particulate  biogenic 
silica  is  produced  through  diatom  mortality.  Particulate  biogenic  silica  under¬ 
goes  dissolution  to  available  silica  or  settles  to  the  bottom  sediments. 

Chemical  Oxygen  Demand 

Chemical  oxygen  demand  is  modeled  as  the  concentration  of  reduced  substances 
that  are  oxidizable  by  inorganic  means.  The  primary  component  of  chemical 
oxygen  demand  is  sulfide  released  from  sediments.  Oxidation  of  sulfide  to 
sulfate  may  remove  substantial  quantities  of  dissolved  oxygen  from  the  water 
column. 

Dissolved  Oxygen 

Dissolved  oxygen  is  required  for  the  existence  of  higher  life  forms.  Oxygen 
availability  determines  the  distribution  of  organisms  and  the  flow  of  energy  and 
nutrients  in  an  ecosystem.  Dissolved  oxygen  is  a  central  component  of  the 
water  quality  model. 

Total  Active  Metal 


Both  phosphate  and  dissolved  silica  sorb  to  inorganic  solids  -  primarily  iron  and 
manganese.  Sorption  and  subsequent  settling  is  one  pathway  for  removal  of 
phosphate  and  silica  from  the  water  column.  The  sediment  model  does  not  in¬ 
clude  iron  and  manganese.  Instead,  total  active  metal  is  defined  as  the  total 
concentration  of  metals  that  are  active  in  phosphate  and  silica  adsorption  and 
transport.  Total  active  metal  is  partitioned  between  paniculate  and  dissolved 
phases  by  an  oxygen-dependent  partition  coefficient. 

Salinity 

Salinity  is  a  conservative  tracer  that  provides  verification  of  the  transport  com¬ 
ponent  of  the  model  and  facilitates  examination  of  conservation  of  mass.  Salini¬ 
ty  also  influences  the  dissolved  oxygen  saturation  concentration  and  is  used  in 
the  determination  of  kinetics  constants  that  differ  in  saline  and  fresh  water. 
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Temperature 


Biochemical  reactions  are  a  function  of  temperature.  Reaction  rates  increase 
with  increasing  temperature.  Extreme  temperatures  cause  algal  mortality. 

Kinetic  Formulations 

The  remainder  of  this  appendix  details  the  kinetics  of  the  mass-conservation 
equation  for  each  state  variable.  Parameters  are  defined  where  they  first  ap¬ 
pear.  For  consistency  with  reported  rate  coefficients,  kinetic  reactions  are 
defined  with  a  temporal  dimension  of  days.  Within  the  CE-QUAL-ICM  code, 
kinetics  sources  and  sinks  rates  are  converted  to  seconds  for  the  mass-conser¬ 
vation  equation. 

Algae 

Algae  play  a  central  role  in  the  carbon  and  nutrient  cycles.  Equations  govern¬ 
ing  the  three  algal  groups  are  similar  with  differences  among  groups  expressed 
through  the  magnitudes  of  coefficients  in  the  equations.  The  letter  "x"  is  used 
as  a  "wild  card"  in  the  equations  and  is  replaced  with  a  letter  that  indicates  a 
specific  algal  group.  Algal  groups  are  indicated  by: 

c  =  cyanobacteria 
d  =  diatoms 
g  =  green  algae 

Sources  and  sinks  of  algae  are: 

Growth  (production) 

Basal  metabolism 

Predation 

Settling 

The  governing  equation  for  algal  biomass  is: 

i-Bx  =  fpx  -  BMx  -  PRx  -  WSx-il  Bx  (A-2) 

6t  [  bz 

where 

Bx  =  algal  biomass,  expressed  as  carbon  (g  C  m^) 
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Px  =  production  (day  ') 

BMx  =  basal  metabolism  (day  ') 

PRx  =  predation  (day  ') 

WSx  =  settling  velocity  (m  day  ') 
z  =  vertical  coordinate  (m) 

Production.  Production  by  phytoplankton  is  determined  by  nutrient  availabilty, 
light  intensity,  and  temperature.  The  effects  are  considered  to  be  multiplica¬ 
tive: 


Px  =  PMx  f(N)  f(I)  f(T)  (^-3) 


where 


PMx  =  production  under  optimal  conditions  (day ') 
f(N)  =  effect  of  suboptimal  nutrient  concentration  (0  _<  f  ^  1) 
f(I)  =  effect  of  suboptimal  illumination  (0  _<  f  _<  1) 
f(T)  =  effect  of  suboptimal  temperature  (0  _<  f  _<  1) 

Cyanobacteria  are  treated  as  freshwater  organisms  that  undergo  rapid  mortality 
in  salt  water.  This  effect  is  included  by  a  salinity  toxicity  term  in  the  cyanobac¬ 
teria  production  equation: 


Pc  =  PMc  f(N)  f(I)  f(T)  f(S) 


(A-4) 


where 


f(S)  —  salinity  effect  on  cyanobacteria  production  (0  ^  1) 

Nutrients.  Carbon,  nitrogen,  and  phosphorus  are  the  primary  nutrients  re¬ 
quired  for  algal  growth.  Diatoms  also  require  silica.  Inorganic  carbon  limita¬ 
tion  is  not  considered  in  the  model.  The  effects  of  the  remaining  nutrients  on 
growth  are  described  by  Monod  kinetics  (Monod,  1949).  In  this  formulation, 
growth  is  dependent  upon  nutrient  availability  at  low  nutrient  concentrations  but 
is  independent  of  nutrients  at  high  concentrations.  A  key  parameter  in  the 
formulation  is  the  half-saturation  concentration.  In  this  formulation,  the  growth 
rate  is  half  the  maximum  when  available  nutrient  concentration  equals  the  half¬ 
saturation  concentration.  Liebig’s  "law  of  the  minimum"  (Odum,  1971)  states 
that  growth  is  determined  by  the  nutrient  in  least  supply.  For  cyanobacteria 
and  greens: 
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f(N)  =  minimum 


(A-5) 


NH^  +  NO3  PO,d 

KHnx  +  NH4  +  NO3  ’  KHpx  +  PO,d 

where 

NH4  =  ammonium  concentration  (g  N  m'^) 

NOS  =  nitrate  concentration  (g  N  m'^) 

KHnx  =  nitrogen  uptake  half-saturation  constant  (g  N  m'^) 

P04d  =  dissolved  phosphate  concentration  (g  P  m"^) 

KHpx  =  phosphorus  uptake  half-saturation  constant  (g  P  m"^) 

Some  cyanobacteria,  notably  the  bloom-forming  genus  Anabaena,  utilize  atmo¬ 
spheric  nitrogen  to  supply  nitrogen  requirements.  In  the  model,  a  non-nitrogen 
fixing  algae  is  considered  the  predominant  cyanobacteria.  The  potential  growth 
of  nitrogen  fixers,  such  as  Anabaena,  can  be  estimated  by  removing  the  nitro¬ 
gen  limitation  on  cyanobacteria. 

Diatoms  require  silica  as  well  as  nitrogen  and  phosphorus  for  growth.  For 
diatoms,  the  nutrient  limitation  is  the  minimum  of  the  limitations  expressed  in 
Equation  A-5  or  the  following: 

f(N)  =  _ — _  (A-6) 

KHs  +  SAd 

where 

SAd  =  dissolved  available  silica  concentration  (g  Si  m'^) 

KHs  =  silica  half-saturation  constant  for  uptake  by  diatoms  (g  Si  m'^) 

Light.  Algal  production  increases  as  a  function  of  light  intensity  until  an  opti¬ 
mal  intensity  is  reached.  Beyond  the  optimal  intensity,  production  declines  as 
intensity  increases.  Steele’s  equation  (DiToro,  et  al.,  1971)  describes  this  phe¬ 
nomenon: 

f(I)  =  J-e‘’'^ 

Is 

where 

I  =  illumination  rate  (Langleys  day  ') 

Is  =  optimal  illumination  (Langleys  day  ') 
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Steele’s  equation  describes  the  instantaneous  light  limitation  at  a  point  in  space. 
However,  the  model  computes  limitation  integrated  over  discrete  time  intervals 
and  aggregated  spatially  by  model  segments.  In  the  model,  Steele’s  equation  is 
integrated  over  one  day  and  is  averaged  over  the  thickness  of  each  model  seg¬ 
ment.  This  interval  does  not  preclude  computation  steps  less  than  a  day  but 
frees  the  model  from  accounting  for  illumination  in  "real  time."  Daily  averag¬ 
ing  precludes  computation  of  diurnal  fluctuations  in  algal  production.  This 
restriction  is  not  severe  since  the  equations  for  algal  growth  are  not  appropriate 
for  short  time  scales. 

Assuming  light  intensity  declines  exponentially  with  depth,  the  integrated, 
averaged  form  of  Steele’s  equation  is: 

f(I)  =  -  ei  (A-8) 

Kess  Az 

ab  =  -  e  ■  ■*  f  A-9) 

FDIs 


at 


lo 

FD  Is 


g  -  Kess  ZD 


(A-10) 


where 


lo  —  daily  illumination  at  water  surface  (Langleys  day’) 

FD  =  fractional  daylength  (0  ^  FD  ^  1) 

Kess  =  total  light  attenuation  coefficient  (m  ') 

Az  =  model  segment  thickness  (m) 

ZD  =  distance  from  water  surface  to  top  of  model  segment  (m) 

Light  attenuation  in  the  water  column  is  composed  of  two  fractions:  a  back¬ 
ground  value  dependent  on  water  color  and  concentration  of  suspended  parti¬ 
cles,  and  extinction  due  to  light  absorption  by  chlorophyll: 

Kess  =  Keb  +  Kechl  Y  _ 1 _ Bx  (A- 11) 

CChlx 


where 


Keb  =  background  light  attenuation  (m  ') 

Kechl  =  chlorophyll  ’a’  light  attenuation  coefficient  (m-  mg  ') 
CChlx  =  algal  group  x  carbon-to-chlorophyll  ratio  (g  C  mg  '  chi) 
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Optimal  illumination  for  photosynthesis  depends  on  algal  taxonomy,  duration  of 
exposure,  temperature,  nutrient  status,  and  previous  acclimation.  Variations  in 
optimal  illumination  are  largely  due  to  adaptations  by  algae  intended  to  maxi¬ 
mize  production  in  a  variable  environment.  Steele  (1962)  noted  the  result  of 
adaptations  is  that  optimal  illumination  is  a  consistent  fraction  50%)  of  daily 
illumination.  Kremer  and  Nixon  (1978)  reported  an  analogous  finding  that 
maximum  algal  production  occurs  at  a  constant  depth  (=  Im)  in  the  water 
column.  Their  approach  is  adopted  here  so  that  optimal  illumination  is  express¬ 
ed: 

Isx  =  loavg  e  ■  (A*  12) 


where 

loavg  =  adjusted  surface  illumination  (Langleys  day  ') 

Doptx  =  maximum  algal  production  depth  (m) 

A  minimum,  Ismin,  is  specified  for  optimal  illumination  so  that  algae  do  not 
thrive  at  extremely  low  light  levels.  The  time  required  for  algae  to  adapt  to 
changes  in  illumination  is  accounted  for  by  computing  Is  based  on  a  time- 
weighted  average  of  daily  illumination: 


loavg  =  0.7  lo  +  0.211  +  0.1 12 


(A-13) 


where 

11  =  daily  illumination  one  day  preceding  model  day  (Langleys  day  ') 

12  =  daily  illumination  two  days  preceding  model  day 

(Langleys  day  ') 

Temperature.  Algal  production  increases  as  a  function  of  temperature  until  an 
optimum  temperature  is  reached.  Above  the  optimum,  production  declines  until 
a  lethal  temperature  is  reached.  Numerous  functional  representations  of  tem¬ 
perature  effects  are  available.  Inspection  of  growth  versus  temperature  curves 
indicates  a  function  similar  to  a  Gaussian  probability  curve  provides  a  good  fit 
to  observations: 


f(T)  =  g-  KTgxI  (T  -  Tmx)'  J  <  Jjjjx 

_  g  -  KTgx2  amx  -  T)"  yy^gn  T  >  Tiux 


(A-14) 
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where 


T  =  temperature  (C°) 

Tmx  =  optimal  temperature  for  algal  growth  (C®) 

KTgxl  =  temperature  effect  below  TMx  on  growth  (C°  ‘) 

KTgx2  =  temperature  effect  above  TMx  on  growth  (C°  -) 

Salinity  Toxicity.  The  effect  of  salinity  on  freshwater  cyanobacteria  is  repre¬ 
sented  by  an  empirical  equation: 


f(S)  = 


Stox^ 
Stox^  + 


(A-15) 


where 


Stox  =  salinity  which  Microcystis  growth  is  halved  (ppt) 

Settling.  In  late  winter  and  early  spring,  diatom  populations  are  assumed  to  be 
characterized  by  large  species  with  high  settling  velocities.  At  the  end  of 
spring,  large  diatoms  are  assumed  to  be  replaced  by  populations  of  smaller 
individuals  with  lower  settling  velocities.  Diatom  settling  is  represented  in  the 
model  as: 

WSd  =  WSdb  +  PrdvalWSds 

where 

WSdb  =  base  diatom  settling  velocity  (m  day  ') 

WSds  =  large  diatoms  enhanced  settling  velocity  (m  day  *) 

Prdval  =  piecewise  function  used  to  specify  seasonal  diatom  settling 
velocity 

Basal  Metabolism.  In  the  model,  basal  metabolism  is  the  sum  of  all  internal 
processes  that  decrease  algal  biomass.  A  portion  of  metabolism  is  respiration 
which  may  be  viewed  as  a  reversal  of  production.  In  respiration,  carbon  and 
nutrients  are  returned  to  the  environment  accompanied  by  the  consumption  of 
dissolved  oxygen.  A  second  internal  sink  of  biomass  is  the  excretion  of  dis¬ 
solved  organic  carbon. 

Respiration  cannot  proceed  in  the  absence  of  oxygen.  Basal  metabolism  cannot 
decrease  in  proportion  to  oxygen  availability  or  algae  would  approach  immortal¬ 
ity  under  anoxic  conditions.  To  solve  this  dilemma,  basal  metabolism  is  con¬ 
sidered  to  be  independent  of  dissolved  oxygen  concentration  but  the  distribution 
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of  metabolism  between  respiration  and  excretion  is  oxygen-dependent.  When 
oxygen  is  freely  available,  respiration  is  a  large  fraction  of  the  total.  When 
oxygen  is  restricted,  excretion  becomes  dominant.  Formulation  of  this  process 
is  detailed  in  the  text  that  describes  algal  effects  on  carbon  and  dissolved  oxy¬ 
gen. 

Basal  metabolism  is  commonly  considered  to  be  an  exponentially  increasing 
function  of  temperature: 


BMx  =  BMrx  e 


(A-17) 


where 

BMrx  =  metabolic  rate  at  Trx  (day‘) 

KTbx  =  temperature  effect  on  metabolism  (C°  ') 

Trx  =  reference  temperature  for  metabolism  (C”) 

Predation.  Detailed  specification  of  predation  rates  requires  predictive  model¬ 
ing  of  zooplankton  biomass  and  activity.  At  present,  zooplankton  are  not 
included  in  the  model.  Consequently,  a  constant  predation  rate  is  specified. 
The  present  formulation  assumes  zooplankton  biomass  is  a  constant  fraction  of 
algal  biomass.  Zooplankton  activity  is  assumed  to  be  influenced  by  temperature 
and  is  taken  into  account  by  incorporating  an'  exponential  temperature  relation¬ 
ship  into  the  predation  term.  The  predation  formulation  is  identical  to  basal 
metabolism.  The  difference  in  predation  and  basal  metabolism  lies  in  the  distri¬ 
bution  of  the  end  products  of  these  processes: 

PRx  =  BPRxe'^^’''^  (A- 18) 


where 

BPRx  =  predation  rate  at  Trx  (day ') 

Effect  of  Algae  on  Organic  Carbon.  During  production  and  respiration,  algae 
primarily  take  up  and  produce  carbon  dioxide,  an  inorganic  form  not  considered 
in  the  model.  A  small  fraction  of  basal  metabolism  is  exuded  as  dissolved 
organic  carbon.  In  the  model,  this  fraction  increases  as  dissolved  oxygen 
becomes  scarce.  Algae  also  produce  organic  carbon  through  the  effects  of 
predation.  Zooplankton  take  up  and  redistribute  algal  carbon  through  grazing, 
assimilation,  respiration,  and  excretion.  Since  zooplankton  are  not  included  in 
the  model,  routing  of  algal  carbon  through  zooplankton  is  simulated  by  empiri- 
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cal  distribution  coefficients.  The  effects  of  algae  on  organic  carbon  are  ex¬ 
pressed  as; 

i-DOC  =  |fCDx  +  (1 -FCDx) _ _ IbMx  +  FCDPPRx  Bx  (A- 19) 

6t  L  [  KHrx  +  DO  J 

—  LPOC  =  FCLP  PRx  Bx  (A-20) 

6t 


1  RPOC  =  FCRP  PRx  Bx  (A-21) 

ot 

where 

DOC  =  dissolved  organic  carbon  concentration  (g  C  m'^) 

DO  =  dissolved  oxygen  concentration  (g  O2  nr^) 

LPOC  =  labile  particulate  organic  carbon  concentration  (g  C  m'^) 

RPOC  =  refractory  particulate  organic  carbon  concentration  (g  C  m'^) 

FCDx  =  fraction  of  basal  metabolism  exuded  as  dissolved  organic  carbon 
KHrx  =  half-saturation  concentration  for  algal  dissolved  organic  carbon 
excretion  (g  O2  m‘^) 

FCDP  =  fraction  of  dissolved  organic  carbon  produced  by  predation 
FCLP  =  fraction  of  labile  particulate  carbon  produced  by  predation 
FCRP  =  fraction  of  refractory  particulate  carbon  produced  by  predation 

The  sum  of  the  three  predation  fractions  must  equal  unity. 

Effect  of  Algae  on  Phosphorus.  Algae  take  up  dissolved  phosphate  during 
production  and  release  dissolved  phosphate  and  organic  phosphorus  through 
mortality.  As  with  carbon,  the  fate  of  algal  phosphorus  released  by  metabolism 
and  predation  is  represented  by  distribution  coefficients.  Since  the  total  phos¬ 
phate  state  variable  includes  both  intra-  and  extracellular  phosphate,  no  explicit 
representation  of  the  effects  of  algae  on  phosphate  is  necessary.  Distribution  of 
total  phosphate  is  determined  by  partition  coefficients  as  detailed  in  the  "PHOS¬ 
PHORUS"  section  of  this  appendix.  The  equations  that  express  the  effects  of 
algae  on  organic  phosphorus  are: 

^  DOP  =  (BMx  FPDx  +  PRx  FPDP)  APC  Bx  (A-22) 

ot 
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A  LPOP  =  (BMx  FPLx  +  PRx  FPLP)  APC  Bx 


(A-23) 


1  RPOP  =  (BMx  FPRx  +  PRx  FPRP)  APC  Bx  (A-24) 


where 


DOP 

LPOP 

RPOP 

APC 

FPDx 

FPLx 

FPRx 

FPDP 

FPLP 

FPRP 


=  dissolved  organic  phosphorus  (g  P  m'^) 

=  labile  particulate  organic  phosphorus  (g  P  m 
=  refractory  particulate  organic  phosphorus  (g  P  m'^) 

==  phosphorus-to-carbon  ratio  of  all  algal  groups  (g  P  g  '  C) 

=  fraction  of  dissolved  organic  phosphorus  produced  by  metabo¬ 
lism 

=  fraction  of  labile  particulate  phosphorus  produced  by  metabo¬ 
lism 

=  fraction  of  refractory  particulate  phosphorus  produced  by  metab¬ 
olism 

=  fraction  of  dissolved  organic  phosphorus  produced  by  predation 
=  fraction  of  labile  particulate  phosphorus  produced  by  predation 
=  fraction  of  refractory  particulate  phosphorus  produced  by  preda¬ 
tion 


The  sums  of  the  metabolism  and  respiration  fractions  must  each  be  less  than  or 
equal  to  unity. 

Effect  of  Algae  on  Nitrogen.  Algae  take  up  ammonium  and  nitrate  during 
production  and  release  ammonium  and  organic  nitrogen  through  mortality. 
Nitrate  is  internally  reduced  to  ammonium  before  biomass  synthesis  occurs 
(Parsons,  et  al.,  1984).  Trace  concentrations  of  ammonium  inhibit  nitrate 
reduction  so  that,  in  the  presence  of  ammonium  and  nitrate,  ammonium  is 
utilized  first.  The  preference  of  algae  for  ammonium  can  be  expressed  empiri¬ 
cally  (Thomann  and  Fitzpatrick,  1982): 


PNx  =  NH, 


NO, 


+  NH 


(KHnx  +  NH,)(KHnx  +  NOj) 
KHnx 


(A-25) 


(NH,  +  NO,)  (KHnx  +  NO,) 
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where 


PNx  =  algal  preference  for  ammonium  uptake  (0  ^  PNx  ^  1) 

The  ammonium  preference  function  has  two  limiting  values.  When  nitrate  is 
absent,  the  preference  for  ammonium  is  unity.  When  ammonium  is  absent,  the 
preference  is  zero.  In  the  presence  of  ammonium  and  nitrate,  the  preference 
depends  on  the  abundance  of  both  forms  relative  to  the  half-saturation  constant 
for  nitrogen  uptake.  When  both  ammonium  and  nitrate  are  abundant,  the 
preference  for  ammonium  approaches  unity.  When  ammonium  is  scarce  but 
nitrate  is  abundant,  the  preference  decreases  in  magnitude  and  a  significant 
fraction  of  algal  nitrogen  comes  from  nitrate. 

Algal  nitrogen  releases  by  metabolism  and  predation  are  represented  by  distri¬ 
bution  coefficients.  The  effects  of  algae  oh  the  nitrogen  state  variables  are 
expressed  as: 

A  NH^  =  (BMx  FNIx  +  PRx  FNIP  -  PNx  Px)  ANCx  Bx  (A-26) 
ot 


i-NOj  =  (PNx  -  l)PxANCxBx  (A-27) 

ot 

DON  =  (BMx  FNDx  +  PRx  FNDP)  ANCx  Bx  (A-28) 

Ot 

LPON  =  (BMx  FNLx  +  PRx  FNLP)  ANCx  Bx  (A-29) 
Ot 

A  RPON  =  (BMx  FNRx  +  PRx  FNRP)  ANCx  Bx  (A-30) 

ot 


where 

DON  =  dissolved  organic  nitrogen  concentration  (g  N  m'^) 

LPON  =  labile  particulate  organic  nitrogen  concentration  (g  N  m 
RPON  =  refractory  particulate  organic  nitrogen  concentration  (g  N  m'^) 
ANCx  =  nitrogen-to-carbon  ratio  of  algae  (g  N  g  ‘  C) 

FNIx  =  fraction  of  inorganic  nitrogen  produced  by  metabolism 
FNDx  =  fraction  of  dissolved  organic  nitrogen  produced  by  metabolism 
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FNLx  =  fraction  of  labile  particulate  nitrogen  produced  by  metabolism 
FNRx  =  fraction  of  refractory  particulate  nitrogen  produced  by  metab¬ 
olism 

FNIP  =  fraction  of  inorganic  nitrogen  produced  by  predation 
FNDP  =  fraction  of  dissolved  organic  nitrogen  produced  by  predation 
FNLP  =  fraction  of  labile  particulate  nitrogen  produced  by  predation 
FNRP  =  fraction  of  refractory  particulate  nitrogen  produced  by  predation 

The  sums  of  the  metabolism  fractions  and  the  predation  fractions  must  each 
equal  unity. 

Algal  Stoichiometry.  Algal  biomass  is  quantified  in  units  of  carbon.  In  order 
to  express  the  effects  of  algae  on  nitrogen  and  phosphorus,  the  ratios  of  nitro- 
gen-to-carbon  and  phosphorus-to-carbon  in  algal  biomass  must  be  specified. 
Global  mean  values  of  these  ratios  are  well  known  (Redfield,  et  al.,  1966). 
Algal  composition  is  known  to  vary  as  a  function  of  nutrient  availability.  As 
nitrogen  and  phosphorus  become  scarce,  algae  adjust  their  composition  so  that 
smaller  quantities  of  these  nutrients  are  required  to  produce  carbonaceous  bio¬ 
mass  (Droop,  1973;  DiToro,  1980;  Parsons  et  al.,  1984).  The  ratio  is  describ¬ 
ed  by  an  empirical  approximation: 


APC  =  - 1 _ 

PCprml  +  PCprm2  e ' 


(A-31) 


where 

PCprml  =  minimum  carbon-to-phosphorus  ratio  (g  C  g  ‘  P) 

PCprm2  =  difference  between  minimum  and  maximum  carbon-to-phos- 
phorus  ratio  (g  C  g  '  P) 

PCprm3  =  effect  of  dissolved  phosphate  concentration  on  carbon-to-phos- 
phorus  ratio  (m^  g  '  P) 

Effect  of  Algae  on  Dissolved  Oxygen.  Algae  produce  oxygen  during  photo¬ 
synthesis  and  consume  oxygen  through  respiration.  The  quantity  produced 
depends  on  the  form  of  nitrogen  utilized  for  growth.  More  oxygen  is  produc¬ 
ed,  per  unit  of  carbon  fixed,  when  nitrate  is  the  algal  nitrogen  source  than  when 
ammonium  is  the  source.  Equations  describing  algal  uptake  of  carbon  and 
nitrogen  and  production  of  dissolved  oxygen  (Morel,  1983)  are: 
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106  CO^  +  16  NH;  +  H,PO;  +  106  H,0  (A-32) 

protoplasm  +  106  0,  +  15  H* 

106  COj  +  16N03"  +  H^PO;  +  122  Hp  +  17  (A-33) 

protoplasm  +  138  0, 


When  ammonium  is  the  nitrogen  source,  one  mole  of  oxygen  is  produced  per 
mole  of  carbon  dioxide  fixed.  When  nitrate  is  the  nitrogen  source,  1.3  moles 
oxygen  are  produced  per  mole  of  carbon  dioxide  fixed. 


The  equation  that  describes  the  effect  of  algae  on  dissolved  oxygen  is: 


i-DO 

6t 


(1.3  -  0.3PNx)Px  - 


DO 


KHrx  +  DO 


BMx 


(A-34) 
AOCR  Bx 


where 

AOCR  =  dissolved  oxygen-to-carbon  ratio  in  respiration  (2.67  g  O2  g  '  C) 

The  magnitude  of  AOCR  is  derived  from  a  simple  representation  of  the  respira¬ 
tion  process: 


CHp  +  O2  =  COj  +  Hp  (A-35) 

The  quantity  (1.3  -  0.3  PNx)  is  the  photosynthesis  ratio  and  expresses  the  molar 
quantity  of  oxygen  produced  per  mole  of  carbon  fixed.  The  photosynthesis 
ratio  approaches  unity  as  the  algal  preference  for  ammonium  approaches  unity. 

Organic  Carbon 

Organic  carbon  undergoes  numerous  transformations  in  the  water  column.  The 
model  carbon  cycle  consists  of  the  following: 

Phytoplankton  production 
Phytoplankton  exudation 
Predation  on  phytoplankton 
Dissolution  of  particulate  carbon 
Heterotrophic  respiration 
Denitrification 
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Settling 


Algal  production  is  the  primary  carbon  source  although  carbon  also  enters  the 
system  through  external  loadings.  Predation  of  algae  releases  particulate  and 
dissolved  organic  carbon  to  the  water  column.  A  fraction  of  the  particulate 
organic  carbon  undergoes  first-order  dissolution  to  dissolved  organic  carbon. 
The  remainder  settles  to  the  sediments.  Dissolved  organic  carbon  produced  by 
phytoplankton  exudation,  predation,  and  dissolution  is  respired  or  denitrified  as 
a  first-order  rate.  No  carbon  is  recycled  from  the  sediments  to  the  water  col¬ 
umn  although  oxygen  demand  created  by  carbon  diagenesis  is  included  in  the 
model. 

Dissolution  and  Respiration  Rates.  Dissolution  and  respiration  rates  depend 
on  the  availability  of  carbonaceous  substrate  and  on  heterotrophic  activity. 
Heterotrophic  activity  and  biomass  have  been  correlated  with  algal  activity  and 
biomass  across  a  wide  range  of  natural  systems  (Bird  and  Kalff,  1984;  Cole,  et 
al.,  1988).  Consequently,  algal  biomass  can  be  incorporated  into  dissolution 
and  respiration  rate  formulations  as  a  surrogate  for  heterotrophic  activity.  The 
correlation  between  algae  and  heterotrophs  occurs  because  algae  produce  labile 
carbon  that  fuel  heterotrophic  activity.  Dissolution  and  respiration  processes  do 
not  require  the  presence  of  algae  and  may  be  fueled  entirely  by  external  carbon 
inputs.  Representation  of  dissolution  and  respiration  in  the  model  allows  speci¬ 
fication  of  algal-dependent  and  algal-independent  rates: 

Kdoc  =  Kdc  +  Kdcalg  Bx  (A-36) 

x-c,d,g 


where 


Kdoc  =  dissolved  organic  carbon  respiration  rate  (day ') 

Kdc  =  minimum  respiration  rate  (day  ') 

Kdcalg  =  constant  that  relates  respiration  to  algal  biomass  (m^  g  '  C  day  ') 


Klpoc  =  Klc  +  Klcalg  ^  Bx  (A-37) 

x=c,d,g 


where 

Klpoc  =  labile  particulate  organic  carbon  dissolution  rate  (day  ') 

Klc  =  minimum  dissolution  rate  (day  ') 

Klcalg  =  constant  that  relates  dissolution  to  algal  biomass  (m^  g  '  C  day  ') 
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(A-38) 


Krpoc  =  Krc  +  Krcaig  Bx 

x-c,d,g 

where 

Krpoc  =  refractory  particulate  organic  carbon  dissolution  rate  (day  ') 

Krc  =  minimum  dissolution  rate  (day  ') 

Krcaig  =  constant  that  relates  dissolution  to  algal  biomass  (m^  g  '  C  day  ') 

An  exponential  function  relates  dissolution  and  respiration  to  temperature. 

Denitrification.  As  oxygen  is  depleted  from  natural  systems,  oxidation  of 
organic  matter  is  effected  by  the  reduction  of  alternate  oxidants  (referred  to  as 
"alternate  electron  acceptors").  The  sequence  in  which  alternate  acceptors  are 
employed  is  determined  by  the  thermodynamics  of  oxidation-reduction  reac¬ 
tions.  The  first  substance  reduced  in  the  absence  of  oxygen  is  nitrate.  A 
representation  of  the  denitrification  reaction  can  be  obtained  by  balancing  stan¬ 
dard  half-cell  redox  reactions  (Stumm  and  Morgan,  1981): 


4  NO3"  +  4  H  *  +5  CHp  2  +  7  Hp  +  5  CO^  (A-39) 

Equation  A-39  describes  the  stoichiometry  of  the  denitrification  reaction.  The 
kinetics  in  the  model  are  first-order.  The  dissolved  organic  carbon  respiration 
rate,  Kdoc,  is  modified  so  that  significant  decay  via  denitrification  occurs  only 
when  nitrate  is  freely  available  and  dissolved  oxygen  is  depleted.  A  parameter 
is  included  so  that  the  anoxic  respiration  rate  is  slower  than  oxic  respiration: 


Denit 


KHodoc  NO3 

KHodoc  +  DO  KHndn  +  NO3 


AANOX  Kdoc 


(A-40) 


where 


Denit  =  dissolved  organic  carbon  denitrification  rate  (day  ') 

AANOX  =  ratio  of  denitrification  to  oxic  carbon  respiration  rate  (0  ^ 
AANOX  ^  1) 

KHodoc  =  dissolved  oxygen  half-saturation  concentration  required  for  oxic 
respiration  (g  O2  m'^) 

KHndn  =  nitrate  half-saturation  concentration  required  for  denitrification 
(g  N  m-^) 


An  exponential  function  relates  denitrification  to  temperature.  Parameter  values 
in  the  function  are  the  same  as  for  dissolved  organic  carbon  respiration. 
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Dissolved  Organic  Carbon.  The  complete  representation  of  all  dissolved 
organic  carbon  sources  and  sinks  in  the  model  is  given  by: 


i-DOC 

=  T 

FCDx+(l -FCDx) 

BMx  +  FCDP  PRx 

at 

x«c,d,g 

KHrx  +  DO 

+  Klpoc  LPOC  +  Krpoc  RPOC  - - - Kdoc  DOC 

^  ^  KHodoc  +  DO 


-Denit  DOC  (A-41) 


Labile  Particnlate  Organic  Carbon.  The  complete  representation  of  all  labile 
particulate  organic  carbon  sources  and  sinks  in  the  model  ecosystem  is  given 
by: 

—  LPOC  =  V  FCLP  PRx  Bx  -  Klpoc  LPOC  -  WSl  —  LPOC 

x-c.d.g  (A-42) 

where 

WSl  =  labile  particulate  settling  velocity  (m  day') 

Refractory  Particulate  Organic  Carbon.  The  complete  representation  of  ail 
refractory  particulate  organic  carbon  sources  and  sinks  in  the  model  ecosystem 
is  given  by: 


i-RPOC 

dt 


T  FCRP  PRx  Bx  -  Krpoc  RPOC  -  WSr  —  RPOC 
x-c.d.g  (A-43) 


where 


WSr  =  refractory  particulate  settling  velocity  (m  day') 

Phosphorus 

The  model  phosphorus  cycle  includes  the  following  processes: 


Algal  production  and  metabolism 
Predation 

Hydrolysis  of  particulate  organic  phosphorus 
Mineralization  of  dissolved  organic  phosphorus 
Settling 

Exchange  with  inorganic  solids 
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External  loads  provide  the  ultimate  source  of  phosphorus  to  the  system.  Dis¬ 
solved  phosphate  is  incorporated  by  algae  during  growth  and  released  as  phos¬ 
phate  and  organic  phosphorus  through  respiration  and  predation.  A  portion  of 
the  particulate  organic  phosphorus  hydrolyzes  to  dissolved  organic  phosphorus. 
The  balance  settles  to  the  sediments.  Dissolved  organic  phosphorus  is  mineral¬ 
ized  to  phosphate.  A  portion  of  the  phosphate  sorbs  to  inorganic  solids  and 
settles  to  the  sediments.  Within  the  sediments,  particulate  phosphorus  is  miner¬ 
alized  and  recycled  to  the  water  column  as  dissolved  phosphate. 

Effects  on  phosphorus  of  algal  production,  metabolism,  and  predation  have 
already  been  detailed.  Descriptions  of  hydrolysis  and  mineralization  and  of  the 
total  phosphate  system  follow. 

Hydrolysis  and  Mineralization.  In  the  model,  hydrolysis  is  defined  as  the 
process  by  which  particulate  organic  substances  are  converted  to  dissolved 
organic  form.  Mineralization  is  defined  as  the  process  by  which  dissolved 
organic  substances  are  converted  to  dissolved  inorganic  form.  Conversion  of 
particulate  organic  phosphorus  to  phosphate  proceeds  through  hydrolysis  and 
mineralization.  Direct  mineralization  of  particulate  organic  phosphorus  does 
not  occur. 

Mineralization  of  organic  phosphorus  is  mediated  by  the  release  of  nucleotidase 
and  phosphatase  enzymes  by  bacteria  (Ammerman  and  Azam,  1985;  Chrost  and 
Overbeck,  1987)  and  algae  (Matavulj  and  Flint,  1987;  Chrost  and  Overbeck, 
1987;  Boni,  et  al.,  1989).  Since  the  algae  themselves  release  the  enzyme  and 
since  bacterial  abundance  is  related  to  algal  biomass,  the  rate  of  organic  phos¬ 
phorus  mineralization  in  the  model  is  related  to  algal  biomass.  A  remarkable 
property  of  the  enzyme  process  is  that  alkaline  phosphatase  activity  is  inversely 
proportional  to  phosphate  concentration  (Chrost  and  Overbeck,  1987;  Boni,  et 
al.,  1989).  Put  in  different  terms,  when  phosphate  is  scarce,  algae  stimulate 
production  of  an  enzyme  that  mineralizes  organic  phosphorus  to  phosphate. 
This  phenomenon  is  simulated  by  relating  mineralization  to  the  algal  phosphorus 
nutrient  limitation.  Mineralization  is  highest  when  algae  are  strongly  phospho¬ 
rus  limited  and  is  least  when  no  limitation  occurs. 

Expressions  for  mineralization  and  hydrolysis  rates  are; 
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Kdop  =  Kdp 


(A-44) 


KHp 

^  KHp  +  PO^d 


Kdpalg  Bx 

x=c,d,g 


where 

Kdop  =  dissolved  organic  phosphorus  mineralization  rate  (day  ') 

Kdp  =  minimum  mineralization  rate  (day"') 

KHp  =  mean  half-saturation  constant  for  algal  phosphorus  uptake  = 
(KHpc  -I-  KHpd  -I-  KHpg)/3  (g  P  m'^) 

Kdpalg  =  constant  that  relates  mineralization  to  algal  biomass  (m^  g '  C 
day') 

where 

Klpop  =  labile  particulate  phosphorus  hydrolysis  rate  (day') 

Kip  =  minimum  hydrolysis  rate  (day ') 

Klpalg  =  constant  that  relates  hydrolysis  to  algal  biomass  (m^  g  ‘  C  day  ') 

Krpop  =  Krp  +  ™  ^  Y 

KHp  +  PO^d 

where 

Krpop  =  refractory  particulate  phosphorus  hydrolysis  rate  (day  ') 

Krp  =  minimum  hydrolysis  rate  (day  ') 

Krpalg  =  constant  that  relates  hydrolysis  to  algal  biomass  (m^  g  ‘  C  day  ') 

An  exponential  function  relates  mineralization  and  hydrolysis  rates  to  tempera¬ 
ture. 


When  the  nutrient  concentration  greatly  exceeds  the  half-saturation  concentra¬ 
tion  for  algal  uptake,  the  rate  roughly  equals  the  minimum.  Algal  biomass  has 
little  influence.  As  phosphorus  becomes  scarce  relative  to  the  half-saturation 
concentration,  the  rate  increases.  The  magnitude  of  increase  depends  on  algal 
biomass. 

The  Total  Phosphate  System.  One  fraction  of  total  phosphorus  in  the  water 
column  is  phosphorus  incorporated  in  algal  biomass.  This  fraction  is  computed 
in  the  model  as  the  product  of  algal  biomass  and  APC,  the  phosphorus-to-car- 
bon  ratio.  In  the  environment,  algae  adjust  their  phosphorus  content  in  re- 
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sponse  to  external  conditions.  Algal  phosphorus  content  is  high  when  external 
phosphorus  is  abundant  and  phosphorus  content  is  low  when  phosphorus  is 
scarce.  The  adaptation  of  algae  to  their  environment  indicates  the  phosphorus- 
to-carbon  ratio  should  be  a  variable  in  the  model.  However,  treatment  of  the 
ratio  as  a  variable  greatly  complicates  computation  of  phosphorus  transport  due 
to  the  mixture  of  algal  masses  of  different  composition.  The  complication  is 
avoided  if  intracellular  and  extracellular  phosphorus  are  treated  and  transported 
as  a  single  state  variable.  Intracellular  and  extracellular  concentrations  are 
determined  by  equilibrium  partitioning  of  their  sum. 

Adsorption  is  the  process  in  which  ions  or  molecules  are  attracted  to  the  surface 
of  a  solid.  Phosphate  ions  exhibit  a  strong  adsorption  to  particulate  species  of 
metals  including  iron  and  manganese.  Particulate  and  dissolved  fractions  are 
determined  by  equilibrium  partitioning  of  their  sum. 

The  model  phosphate  state  variable  is  defined  as  the  sum  of  dissolved  phos¬ 
phate,  sorbed  phosphate,  and  algal  phosphorus  content  as: 

PO^t  =  PO,d  +  PO^  +  PO^a  (A-47) 


where 

P04t  =  total  phosphate  (g  P  m'^) 

P04d  =  dissolved  phosphate  (g  P  m 

PO4P  =  particulate  (sorbed)  phosphate  (g  P  m 

P04a  =  algal  phosphorus  (g  P  m‘^) 

Exchange  With  Particulate  Metals.  Detailed  treatment  of  iron  and  manganese 
is  not  included  in  the  model.  Instead,  the  state  variable  total  active  metal  is 
defined  as  the  sum  of  all  metals  that  act  as  sorption  sites.  Phosphate  sorbs  only 
to  the  particulate  fraction  of  the  total  metal.  Therefore,  the  total  metal  is  parti¬ 
tioned  into  particulate  and  dissolved  fractions  via  an  equilibrium  partition  coef¬ 
ficient.  Extracellular  phosphate  is  partitioned  by  a  linear  sorption  isotherm  so 
that  dissolved  and  particulate  fractions  are  computed  as: 
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(A-48) 


PO,d 


1 

1  +  Kadpo4TAMp 


(PO,t  -  PO,a) 


1  +  Kadpo4TAMp 


where 

Kadpo4  =  sorbed  vs.  dissolved  phosphate  partition  coefficient  (m^  mof) 
TAMp  =  particulate  total  active  metal  (mol  m'^) 

Computation  of  Algal  Phosphorus.  Algal  phosphorus  is  defined  as; 

PO,a  =  APC  Bx  (A-50) 

x-c,d,g 


The  phosphorus-to-carbon  ratio  is  calculated  by  the  empirical  function  expressed 
in  Equation  A-3 1 . 

The  equations  4-31,  A-48,  A-49,  A-50,  form  an  implicit  system.  The  implicit 
nature  of  the  system  can  be  seen  by  substituting  the  expressions  for  dissolved 
phosphate  (Equation  A-48)  and  algal  phosphate  (Equation  A-50)  into  the  expres¬ 
sion  for  phosphorus-to-carbon  ratio  (Equation  A-31); 

APC  =  - 1 - -  (A-51) 

PCprml  +  PCprm2  e '  ^ 

The  phosphorus-to-carbon  ratio  appears  on  both  sides  of  the  equation  and  can¬ 
not  be  obtained  explicitly.  (For  simplicity,  Kadpo4  =  0  in  this  example.) 
Within  the  model  code,  the  phosphorus-to-carbon  ratio  is  obtained  by  iterative 
solution  of  Equation  A-51  in  each  control  volume.  From  the  phosphorus-to- 
carbon  ratio,  algal  phosphorus  is  computed  (Equation  A-50).  Next,  extra¬ 
cellular  phosphate  is  partitioned  into  dissolved  and  particulate  fractions  (Eqs.  A- 
48  and  A-49). 

Phosphate.  Once  the  interactions  of  dissolved,  particulate,  and  algal  phosphate 
are  made  explicit,  the  balance  of  the  equations  describing  phosphorus  are 
straightforward  summations  of  previously-described  sources  and  sinks: 
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(A-52) 


1  PO,t  =  -  y  WSx  —  APC  Bx  -  WSs  —  POi)  +  Kdop  DOP 
where 

WSs  =  particulate  metal  settling  velocity  (m  day  ') 

Algal  uptake  and  release  of  phosphate  represents  an  exchange  of  phosphate 
fractions  rather  than  a  phosphate  source  or  sink.  Consequently,  no  algal  source 
or  sink  terms  are  included  in  the  phosphate  mass-conservation  equation.  The 
settling  terms  are  required  to  represent  the  settling  of  particulate  phosphate 
incorporated  in  algal  biomass  or  sorbed  to  particles. 

Dissolved  Organic  Phosphorus. 

-  DOP  =  y  (BMx  FPDx  +  PRx  FPDP)  APC  Bx 
6t  (A-53) 

+  Klpop  LPOP  +  Krpop  RPOP  -  Kdop  DOP 


Labile  Particulate  Organic  Phosphorus. 


6t 


LPOP  =  y  (BMx  FPLx  +  PRx  FPLP)  APC  Bx 


X«C,d,g 


-  Klpop  LPOP  -  WSl  _  LPOP 

8z 


(A-54) 


Refractory  Particulate  Organic  Phosphorus. 


6t 


RPOP  =  y  (BMx  FPRx  +  PRx  FPRP)  APC  Bx 


x«c.d,g 


(A-55) 


-  Krpop  RPOP  -  WSr  _  RPOP 

8z 
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Nitrogen 


The  model  nitrogen  cycle  includes  the  following  processes: 

Algal  production  and  metabolism 
Predation 

Hydrolysis  of  particulate  organic  nitrogen 

Mineralization  of  dissolved  organic  nitrogen 

Settling 

Nitrification 

Denitrification 

External  loads  provide  the  ultimate  source  of  nitrogen  to  the  system.  Inorganic 
nitrogen  is  incorporated  by  algae  during  growth  and  released  as  ammonium  and 
organic  nitrogen  through  respiration  and  predation.  A  portion  of  the  particulate 
organic  nitrogen  hydrolyzes  to  dissolved  organic  nitrogen.  The  balance  settles 
to  the  sediments.  Dissolved  organic  nitrogen  is  mineralized  to  ammonium.  In 
an  oxygenated  water  column,  a  fraction  of  the  ammonium  is  subsequently 
oxidized  to  nitrate.  In  anoxic  water,  nitrate  is  lost  to  nitrogen  gas  through 
denitrification.  Particulate  nitrogen  that  settles  to  the  sediments  is  mineralized 
and  recycled  to  the  water  column,  primarily  as  ammonium.  Nitrate  moves  in 
both  directions  across  the  sediment-water  interface,  depending  on  relative  con¬ 
centrations  in  the  water  column  and  sediment  interstices. 

Effects  on  nitrogen  of  algal  production,  metabolism,  and  predation  have  already 
been  detailed.  Descriptions  of  hydrolysis,  mineralization,  nitrification  and 
denitrification  follow. 

Hydrolysis  and  Mineralization.  In  the  model,  particulate  organic  nitrogen  is 
converted  to  dissolved  organic  nitrogen  via  hydrolysis.  Dissolved  organic 
nitrogen  is  converted  to  ammonium  through  mineralization.  Conversion  of 
particulate  nitrogen  to  ammonium  proceeds  through  hydrolysis  and  mineraliza¬ 
tion.  Direct  mineralization  of  particulate  nitrogen  does  not  occur.  The  argu¬ 
ment  for  accelerated  hydrolysis  and  mineralization  during  nutrient-limited 
conditions  is  not  as  clear  for  nitrogen  as  for  phosphorus.  However,  the  same 
formulations  are  made  available  for  nitrogen  as  for  phosphorus.  Accelerated 
processes  can  be  activated  or  deactivated  through  parameter  selection.  The 
nitrogen  hydrolysis  and  mineralization  formulations  are: 
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Kdon  =  Kdn  + 


KHn 

KHn  +  NH,  +  NO3 


Kdnalg  ^  Bx 

x=c,d,g 


(A-56) 


where 

Kdon 

Kdn 

KHn 

Kdnalg 


where 

Klpon 

Kin 

Klnalg 


where 

Krpon 

Krn 

Krnalg 


=  dissolved  organic  nitrogen  mineralization  rate  (day  ‘) 

=  minimum  mineralization  rate  (day  ') 

=  algal  nitrogen  uptake  mean  half-saturation  constant  =  (KHnc  + 
KHnd  -I-  KHng)/3  (g  N  m 

=  constant  that  relates  mineralization  to  algal  biomass  (m^  g '  C 
day') 


Klpon  =  Kin  + 


KHn 

KHn  +  NH^  +  NO3 


Klnalg  ^  Bx 

x-c,d.g 


(A-57) 


=  labile  particulate  nitrogen  hydrolysis  rate  (day  ') 

=  minimum  hydrolysis  rate  (day  ') 

=  constant  that  relates  hydrolysis  to  algal  biomass  (m^  g  '  C  day  ') 


Krpon  =  Krn  + 


KHn 

KHn  +  NH,  +  NO3 


Krnalg  Bx 

x-c,d,g 


(A-58) 


=  refractory  particulate  nitrogen  hydrolysis  rate  (day') 

=  minimum  hydrolysis  rate  (day') 

=  constant  that  relates  hydrolysis  to  algal  biomass  (m^  g  '  C  day') 


An  exponential  function  relates  mineralization  and  hydrolysis  rates  to  tempera¬ 
ture. 


Nitriflcation.  Nitrification  is  a  process  mediated  by  specialized  groups  of 
autotrophic  bacteria  that  obtain  energy  through  the  oxidation  of  ammonium  to 
nitrite  and  oxidation  of  nitrite  to  nitrate.  A  simplified  expression  for  complete 
nitrification  (Tchobanoglous  and  Schroeder,  1987)  is; 


NH;  +  2  O2  N03'  +  H3O  +  2  H  *  (A-59) 

The  equation  indicates  that  two  moles  of  oxygen  are  required  to  nitrily  one 
mole  of  ammonium  into  nitrate.  However,  the  simplified  equation  is  not  strict¬ 
ly  true.  Cell  synthesis  by  nitrifying  bacteria  is  accomplished  by  the  fixation  of 
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carbon  dioxide  so  that  less  than  two  moles  of  oxygen  are  consumed  per  mole 
ammonium  utilized  (Wezernak  and  Gannon,  1968). 


The  kinetics  of  complete  nitrification  are  modeled  as  a  function  of  available 
ammonium,  dissolved  oxygen,  and  temperature: 


NT 


DO  NH, 

KHont  +  DO  KHnnt  -+  NH, 

4 


f(T)  NTm 


(A-60) 


where 

NT  =  nitrification  rate  (g  N  m’^  day') 

KHont  =  dissolved  oxygen  half-saturation  constant  required  for  nitrifi¬ 
cation  (g  O2  m’^) 

KHnnt  =  NH4  half-saturation  constant  required  for  nitrification  (g  N  m'^) 
NTm  =  maximum  nitrification  rate  at  optimal  temperature 
(g  N  m'^  day') 

The  kinetics  formulation  incorporates  the  products  of  two  Monod  functions. 
The  first  function  diminishes  nitrification  at  low  dissolved  oxygen  concentra¬ 
tion.  The  second  function  expresses  the  influence  of  ammonium  concentration 
on  nitrification.  When  ammonium  concentration  is  low,  relative  to  KHnnt, 
nitrification  is  proportional  to  ammonium  concentration.  For  NH4  <  < 
KHnnt,  the  reaction  is  approximately  first-order.  (The  first-order  decay  con¬ 
stant  «  NTm/KHnnt.)  When  ammonium  concentration  is  large  relative  to 
KHnnt,  nitrification  approaches  a  maximum  rate.  This  formulation  is  based  on 
a  concept  proposed  by  Tuffey,  et  al.  (1974).  Nitrifying  bacteria  adhere  to 
benthic  or  suspended  sediments.  When  ammonium  is  scarce,  vacant  surfaces 
suitable  for  nitrifying  bacteria  exist.  As  ammonium  concentration  increases, 
bacterial  biomass  increases,  vacant  surfaces  are  occupied,  and  the  rate  of  nitrifi¬ 
cation  increases.  The  bacterial  population  attains  maximum  density  when  all 
surfaces  suitable  for  bacteria  are  occupied.  At  this  point,  nitrification  proceeds 
at  a  maximum  rate  independent  of  additional  increase  in  ammonium  concentra¬ 
tion. 

The  optimal  temperature  for  nitrification  may  be  less  than  peak  temperatures 
that  occur  in  coastal  waters.  To  allow  for  a  decrease  in  nitrification  at  superop- 
timal  temperature,  the  effect  of  temperature  on  nitrification  is  modeled  in  the 
Gaussian  form  of  Equation  A- 14. 
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Effect  of  Nitrification  on  Anunonium. 


A  NH,  =  -  NT  (A-61) 

6t  ' 


Effect  of  Nitrification  on  Nitrate. 

—  NO,  =  NT  (A-62) 

5t  ' 


Effect  of  Nitrification  on  Dissolved  Oxygen. 

DO  =  -  AONT  NT  (A-63) 

6t 


where 

AONT  =  mass  dissolved  oxygen  consumed  per  mass  ammonium-nitrogen 
nitrified  (4.33  g  O2  g'‘  N) 

Effect  of  Denitrification  on  Nitrate.  The  effect  of  denitrification  on  dissolved 
organic  carbon  has  been  described.  Denitrification  removes  nitrate  from  the 
system  in  stoichiometric  proportion  to  carbon  removal  as  determined  by  Equa¬ 
tion  A-39: 

I  NO3  =  -  ANDC  Denit  DOC  (A-64) 


where 

ANDC  =  mass  nitrate-nitrogen  reduced  per  mass  dissolved  organic  carbon 
oxidized  (0.933  g  N  g  '  °C) 

Nitrogen  Mass  Balance  Equations.  The  mass-balance  equations  for  nitrogen 
state  variables  are  written  by  summing  all  previously-described  sources  and 
sinks: 
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Ammonium. 


_  NH,  =  (BMx  FNIx  +  PRx  FNIP  -  PNxPx)  ANCx  Bx  (A-65) 

x-c,d.g 

+  KdonDON  -  NT 


Dissolved  Organic  Nitrogen. 

—  DON  =  Y  FNDx  +  PRx  FNDP)  ANCx  Bx  (A-66) 

K-C,d,g 

+  KlponLPON  +  KrponRPON  -  KdonDON 


Labile  Particulate  Organic  Nitrogen. 


LPON  =  Y  (BMx  FNLx  +  PRx  FNLP)  ANCx  Bx 


Klpon  LPON  -  WSl  _  LPON 
8z 


Refractory  Particulate  Organic  Nitrogen. 

A  RPON  =  Y  (BMx  FNRx  +  PRx  FNRP)  ANCx  Bx 

Ot  x-c  H  p 


Krpon  RPON  -  WSr  _  RPON 
8z 


(A-67) 


(A-68) 


Nitrate. 


—  NO3  =  Y  (BNx  -  l)PxANCxBx 


+  NT  -  ANDC  Denit  DOC 


(A-69) 


Silica 


The  model  incorporates  two  siliceous  state  variables,  available  silica  and  partic¬ 
ulate  biogenic  silica.  For  practical  purposes,  available  silica  is  equivalent  to 
dissolved  silica  although  sorption  of  available  silica  to  inorganic  solids  occurs. 
The  silica  cycle  is  a  simple  one  in  which  diatoms  take  up  available  silica  and 
recycle  available  and  particulate  biogenic  silica  through  the  actions  of  metabo- 
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lism  and  predation.  Particulate  silica  dissolves  in  the  water  column  or  settles  to 
the  bottom.  A  portion  of  the  settled  particulate  biogenic  dissolves  within  the 
sediments  and  returns  to  the  water  column  as  available  silica.  Sources  and 
sinks  represented  are: 


Diatom  production  and  metabolism 
Predation 

Dissolution  of  particulate  to  dissolved  silica 
Settling 

Exchange  with  inorganic  solids 


Available  Silica.  Dissolved  silica  sorb  to  newly-formed  metal  particles  and 
rapidly  settled  to  the  bottom.  To  represent  this  process,  available  silica  is 
partitioned  into  dissolved  and  particulate  fractions  according  to  a  linear  sorption 
isotherm: 


SAd 


1 

1  +  KadsaTAMp 


SA 


(A-70) 


SAp  =  KadsaTAMp 
1  +  KadsaTAMp 


(A-71) 


where 


SA  =  available  silica  concentration  (g  Si  m'^’ 

SAd  =  dissolved  available  silica  (g  Si  m^) 

SAp  =  particulate  available  silica  (g  Si  m 
Kadsa  =  partition  coefficient  of  sorbed  vs.  dissolved  available  silica  (m^ 
mol ') 


The  kinetic  equation  for  available  silica  is: 

4  SA  =  (FSAP  PRd  -  Pd)  ASCd  Bd  -  WSs  —  SAp  /  *  79^ 
ot  Sz 

+  KsuaSU 


where 

SU  =  particulate  biogenic  silica  concentration  (g  Si  m'^) 

FSAP  =  fraction  of  diatom  silica  made  available  by  predation  (0  j< 
FSAP  ^  1) 

ASCd  =  silica-to-carbon  ratio  of  diatoms  (g  Si  g  '  C) 

Ksua  =  particulate  silica  dissolution  rate  (day  ') 
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Particulate  Biogenic  Silica.  The  kinetic  equation  for  particulate  biogenic  silica 
is: 

i-SU  =  (BMd  +  (1  -  FSAP)  PRd)  ASCd  Bd  -  WSs SU  (A-'73) 
6t  6z 

-  KSUA  SU 

An  exponential  function  describes  the  effect  of  temperature  on  silica  dissolution. 


Chemical  Oxygen  Demand 

In  the  model,  chemical  oxygen  demand  is  the  concentration  of  reduced  sub¬ 
stances  that  are  oxidizable  through  inorganic  means.  The  source  of  chemical 
oxygen  demand  in  saline  water  is  sulfide  released  from  sediments.  A  cycle 
occurs  in  which  sulfate  is  reduced  to  sulfide  in  the  sediments  and  reoxidized  to 
sulfate  in  the  water  column.  In  freshwater,  methane  is  released  to  the  water 
column  by  the  sediment  model.  Both  sulfide  and  methane  are  quantified  in 
units  of  oxygen  demand  and  are  treated  with  the  same  kinetic  formulation: 


i-COD  =  - 
6t 


DO 

KHocod  +  DO 


Kcod  COD 


(A-74) 


where 

COD  =  chemical  oxygen  demand  concentration  (g  02-equivalents  m‘^) 
KHocod  =  half-saturation  concentration  of  dissolved  oxygen  required  for 
exertion  of  chemical  oxygen  demand  (g  O2  m'^) 

Kcod  =  oxidation  rate  of  chemical  oxygen  demand  (day*) 

An  exponential  function  describes  the  effect  of  temperature  on  exertion  of 
chemical  oxygen  demand. 


Dissolved  Oxygen 

Sources  and  sinks  of  dissolved  oxygen  in  the  water  column  include: 

Algal  photosynthesis 
Atmospheric  reaeration 
Algal  respiration 
Heterotrophic  respiration 
Nitrification 

Chemical  oxygen  demand 
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Reaeration.  The  rate  of  reaeration  is  proportional  to  the  dissolved  oxygen 
deficit  in  model  segments  that  form  the  air- water  interface: 

A  DO  =  ^  (DOS  -  DO)  (A-75) 

6t  Az 


where 


Kr  =  reaeration  coefficient  (m  day  ') 

DOs  =  dissolved  oxygen  saturation  concentration  (g  O2  m  ^) 

The  surface  renewal  concept,  attributed  to  Danckwerts  by  O’Connor  and  Dob¬ 
bins  (1958),  indicates: 

Kr  =  v/DIR  (^-'76) 

where 

D1  =  molecular  diffusivity  of  oxygen  in  water  (~  1.7  x  10-4  m^  day  ^) 
R  =  surface  renewal  rate  (day  ') 

Specification  of  the  surface  renewal  rate  is  the  fundamental  problem  in  reaera¬ 
tion  theory.  O  Connor  and  Dobbins  (1958)  state  that,  in  isotropic  turbulence, 
surface  renewal  can  be  approximated  as  the  ratio  of  stream  velocity  to  depth. 
The  renewal  rate  is  also  influenced  by  wind  (O’Connor,  1983).  Influences  of 
temperature  (ASCE,  1961)  and  salinity  (Wen,  et  al.,  1984)  on  aeration,  most 
likely  effected  through  changes  in  diffusivity,  have  been  measured.  No  single 
theory  that  unites  all  these  factors  into  a  formulation  of  reaeration  in  an  estuary 
is  available.  The  surface  renewal  concept  is  retained  with  the  renewal  rate 
treated  as  a  calibration  parameter. 


Saturation  dissolved  oxygen  concentrations  decrease  as  temperature  and  salinity 
increase.  An  empirical  formula  that  describes  these  effects  (Genet,  et  al.,  1974) 
is: 


DOs  =  14.5532  -  0.38217  T  +  0.0054258  T^ 


(A-77) 


-  (1.665x  10"'  -  5.866  X  10'*  T  +  9.796  x  10'*  T^)  CL 


where 

CL  -  chloride  concentration  (=  salinity/ 1.80655) 
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Summary  of  Dissolved  Oxygen  Sources  and  Sinks.  The  complete  kinetics  for 
dissolved  oxygen  are; 


Ido  =  T 

x-Tig 


(1.3  -  0.3PNx)Px  - 
DO 


DO 


BMx 


-  AONTNT  - 


KHodoc  +  DO 


KHrx  +  DO 
AOCR  Kdoc  DOC 


AOCR  Bx 
(A-78) 


- £2 _ Kcod  COD  +  iSl  (DOs  -  DO) 

KHocod  +  DO  Az 


Total  Active  Metal 


The  total  active  metal  state  variable  is  the  sum  of  iron  and  manganese  concen¬ 
trations.  Iron  and  manganese  exist  in  particulate  and  dissolved  forms  in  estuar¬ 
ies.  The  prime  determinant  of  the  species  is  dissolved  oxygen.  In  the  oxygen¬ 
ated  water,  total  iron  and  manganese  are  almost  completely  particulate.  Under 
anoxic  conditions,  large  fractions  of  total  iron  and  manganese  are  dissolved 
although  solid-phase  sulfides  and  carbonates  exist  and  may  predominate.  The 
simplified  partitioning  of  particulate  and  dissolved  phases  of  total  active  metal 
employed  here  notes  that  total  active  metal  concentration  must  achieve  a  mini¬ 
mum  level  before  precipitation  occurs.  The  minimum  level  is  a  function  of 
dissolved  oxygen: 

TAMd  =  minimum  (TAMdmx  e"  Kdoiamoo^  TAM)  (A-79) 


TAMp  =  TAM  -  TAMd  (^-80) 


where 


TAM 

TAMd 

TAMp 

TAMdmx 

Kdotam 


total  active  metal  concentration  (mol  m 
dissolved  total  active  metal  (mol  m^) 
particulate  total  active  metal  (mol  m 

solubility  of  total  active  metal  under  anoxic  conditions  (mol  m 
constant  that  relates  total  active  metal  solubility  to  dissolved 
oxygen  concentration  (m^  g '  O2) 


The  origin  of  total  active  metal  in  the  model  is  benthic  sediments.  Since  metal 
release  is  not  explicit  in  the  sediment  model,  release  is  incorporated  into  the 
kinetics  portion  of  the  water-column  model.  Release  is  treated  as  a  spatially- 
uniform,  empirical  function  of  temperature  and  dissolved  oxygen  concentration. 
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The  only  other  source  or  sink  in  the  water  column  is  settling  of  the  particulate 
fraction.  In  the  mass  balance  equation,  the  benthic  source  operates  only  in 
model  cells  that  adjoin  the  bottom; 

i-TAM  =  -  WSs-lxAMp  +  _  BENT  AM 

5t  bz  ^  KHbmf  +  DO  Az 


where 

BENTAM  =  anoxic  total  active  metal  release  rate  release  rate  (mol  m^  day  ') 
KHbmf  =  dissolved  oxygen  concentration  at  which  total  active  metal  re¬ 
lease  is  half  the  anoxic  rate  (g  O2  m^) 

Release  of  metal  from  the  bottom  is  maximum  when  dissolved  oxgyen  is  absent 
from  the  overlying  water.  Release  declines  as  oxygen  concentration  increases 
and  is  negligible  when  dissolved  oxygen  >  >  KHbmf. 

Sediment  release  of  total  active  metal  is  treated  as  an  exponential  function  of 
temperature. 

Salinity 

No  internal  sources  or  sinks  of  salinity  exist. 

T  emperature 

The  only  source  or  sink  of  internal  energy  considered  is  exchange  with  the 
atmosphere.  Although  solar  radiation  can  penetrate  several  meters  into  the 
water  column,  radiation-induced  increases  in  internal  energy  are  assigned  entire¬ 
ly  to  the  surface  model  layer. 

The  internal-energy  equation  can  be  written  as  a  conservation  of  temperature 
equation.  Change  in  temperature  due  to  atmospheric  exchange  is  considered 
proportional  to  the  temperature  difference  between  the  water  surface  and  a 
theoretical  equilibrium  temperature  (Edinger,  et  al.,  1965): 


T  =  . . - . (Te  -  T)  (A-82) 

5t  pCpAz  ’ 


where 
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Te  =  equilibrium  temperature  (°C) 

KT  =  Heat  exchange  coefficient  (W  m'^  °C  ') 

Cp  =  specific  heat  of  water  (4200  W  sec  kg  '  °C"') 
p  =  density  of  water  (1000  kg  m  ^) 
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